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Abstract 

An invariant difFerential-geometric approach to the integrability of (2 + l)-dimensional 
systems of hydrodynamic type, 

ut + Aiu)u^ + B{u)uy = 0, 

is developed. We prove that the existence of special solutions known as 'double waves' is 
equivalent to the diagonalizability of an arbitrary matrix of the two-parameter family 

{kE + A)-\IE + B). 

Since the diagonalizability can be effectively verified by calculating the Haantjes tensor, this 
provides a simple necessary condition for integrability. 
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1 Introduction 



Over the last 20 years there has been a considerable progress in the theory of one-dimensional 
systems of hydrodynamic type + t;(u)u^ = or, in components, 

ul + Vj{u)ui = 0, i,j = l,...,m, (1) 

(the standard summation convention over repeated indices is adopted). Such systems naturally 
occur in applications in gas dynamics, fluid mechanics, chemical kinetics, Whitham averaging 
procedure, differential geometry and topological field theory. We refer to |H1 EZi ISSj -36, ^ 
for a further discussion and references. It has been observed that many particularly important 
examples are diagonalizable, that is, reducible to the Riemann invariant form 

Rl + v\R)Rl = (2) 

where the characteristic speeds f*(R) satisfy the so-called semi-Hamiltonian property (also 
known as the 'richness' condition ^38 ), 



dk = d/dEl^, i ^ j ^ k. We emphasize that the semi-Hamiltonian property © is usually 
automatically satisfied for diagonalizable systems of the 'physical' origin. For instance, a con- 
servative diagonalizable system is necessarily semi-Hamiltonian. We recall that an n-component 
system is said to be conservative if it possesses n conservation laws of hydrodynamic type 
whose densities are functionally independent. It turns out that the additional requirement of the 
diagonalizability implies the existence of an infinity of conservation laws and, hence, the semi- 
Hamiltonian property: see jSH]) or the Appendix to ^H] for a simpler proof. Semi-Hamiltonian 
systems possess infinitely many conservation laws and commuting flows of hydrodynamic type 
and can be linearized by the generalized hodograph method jlUI- Their analytic, differential- 
geometric and Hamiltonian aspects are well-understood by now. 

Remarkably, there exist the effective tensor criteria to verify the diagonalizability and the 
semi-Hamiltonian property without the actual computation of eigenvalues and eigenvectors of 
the matrix w*. Let us first calculate the Nijenhuis tensor of the matrix v*, 

N], = vpu.vl - vld^v] - vl{d^,vl - d^.v^), (4) 

and introduce the Haantjes tensor 

H^, = Niy^vl - N^vlvl - N^.v^j + N^.vlv;. (5) 

For strictly hyperbolic systems the condition of diagonalizablity is given by the following theorem 
(which was stated in jJH] in purely geometric terms as a condition of diagonalizability of a (1,1) 
tensor field). 

Theorem 1 >23^ A hydrodynamic type system 0) with mutually distinct characteristic speeds 
is diagonalizable if and only if the corresponding Haantjes tensor (Ejj is identically zero. 

Since components of the Haantjes tensor can be obtained using any computer algebra package, 
this provides an effective diagonalizability criterion. This criterion has been successfully imple- 
mented in to classify isotherms of adsorption for which the equations of chromatography 
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possess Riemann invariants (since the equations of chromatography are conservative, the semi- 
Hamiltonian property will be automatically satisfied). The same criterion was applied in to 
the Whitham equations governing slow modulations of traveling waves for the generalized KdV 
equation ut + f{u)ux + Uxxx = 0. It was demonstrated that the Whitham equations are diago- 
nalizable (and hence semi-Hamiltonian due to their conservative nature) if and only if = 0. 
Although, for conservative systems, the diagonalizability implies the semi-Hamiltonian property, 
this is not true in general. The tensor object responsible for the semi-Hamiltonian property was 
introduced in j^^j (see the Appendix). 

The present paper aims at the discussion of the extent to which the one-dimensional theory 
carries over to (2 + l)-dimensional quasilinear systems 

Ut + A{u)ux + B{u)uy = 0; (6) 

here u is an m-component column vector and ^(u), B{u) are m x m matrices. Systems of this 
type describe many physical phenomena. In particular, important examples occur in gas dynam- 
ics, shallow water theory, combustion theory, general relativity, nonlinear elasticity, magneto- 
fluid dynamics, etc |27| l6]. Particularly interesting 'integrable' systems of the form Q arise as 
dispersionless limits of multi-dimensional soliton equations j43j . within the method of Whitham 
averaging applied to 'integrable' two-dimensional models |24| I25j , and the R-matrix approach 

The first natural restriction to impose is that all systems arising as one-dimensional limits 
of ©, in particular, the systems + A{u)ux = and ut + B{u)uy = 0, are diagonalizable 
(but not simultaneously, as the matrices A and B do not commute in general). Furthermore, 
applying to (jH)) an arbitrary linear change of the independent variables, 

i = aiit + ai2X + aisy, x = a2it + a22X + a23y , y = a^it + a32X + a^sy , 

we arrive at the transformed system 

u^- + i(u)u£ + ^(u)u^ = 

where 

A = {auE + auA + aisB)-'^{a2iE + 022^ + 023^), 
B = (aiiE + ai2A + arsBy^asiE + 032^ + 033^), 

and E is the m x m identity matrix. Since we want our approach to be invariant under linear 
changes of variables, we require that all matrices of the multi-parameter family 

{aE + bA + cB)-^{aE + bA + cB) (7) 

are diagonalizable and semi-Hamiltonian. These considerations motivate the following defini- 
tions: 

Definition \ A {2 + 1)- dimensional system is said to he diagonalizable if an arbitrary matrix 
of the family Q) is diagonalizable. 

Definition 2 A diagonalizable (2 + 1)- dimensional system is said to be semi-Hamiltonian if an 
arbitrary matrix of the family Q) is semi-Hamiltonian. 

Remarks. One can show that some parameters in are, in fact, redundant: it is sufficient 
to verify the diagonalizability and the semi-Hamiltonian property for an arbitrary matrix in the 
smaller family 

{kE + A)-^{IE + B) (8) 



3 



where k and / are arbitrary constants. Indeed, one can simplify the general matrix Q using 
the fact that the inversion and the addition of a multiple of the identity do not effect the 
diagonalizability. 

We point out that for many systems © the diagonalizability is already sufficiently restrictive 
and implies the semi-Hamitonian property. This is the case, for instance, if the original two- 
dimensional system © is conservative. Indeed, all one-dimensional limits of a multi-dimensional 
conservative system inherit the conservative form, and in (H-l) dimensions the diagonalizability 
is known to imply the semi-Hamiltonian property. 

Finally, we point out that both definitions generalize to multi-dimensional setting (3-1-1, etc) 
in an obvious way. 

Examples of diagonalizable semi-Hamiltonian three-component systems © are discussed in 
Sect. 2. The classification of a special class of diagonalizable three-component conservative 
systems is obtained in Sect. 4. 

An alternative approach to the integrability of multi-dimensional systems of hydrodynamic 
type is based on the method of hydrodynamic reductions. The key element of this construction 
are exact solutions of the system © of the form u(R) = u(i2^, i?") where the Riemann 
invariants R^, solve a pair of commuting diagonal systems 

= A^(R) i?;, = ^^\^) K; (9) 

notice that the number of Riemann invariants is allowed to be arbitrary. Thus, the original 
(2-|-l)-dimensional system (0 is decoupled into a compatible pair of diagonal (l-l-l)-dimensional 
systems. Solutions of this type are known as multiple waves, or nonlinear interactions of n planar 
simple waves (also called solutions with a degenerate hodograph: simple waves (double waves) 
when the number of Riemann invariants equals one (two)). These solutions were extensively 
investigated in gas dynamics and magnetohydrodynamics in a series of publications [HHl [2l IHl IS 
EH ESI 13 1^- Later, they reappeared in the context of the dispersionless KP and Toda hierarchies 
[lailHlIlSlEaEHllMlliaiHl, the theory of integrable hydrodynamic chains [Hll ESI EOl and 
the Laplacian growth problems [2^] • In QS] > it was suggested to call a multi-dimensional system 
integrable if, for arbitrary n, it possesses infinitely many n-component reductions of the form Q 
parametrized by n arbitrary functions of a single argument. It was shown that this requirement 
provides an effective classification criterion. Partial classification results were obtained in |13| 
I14| II 5j . It was demonstrated in |1H II 5j that the method of hydrodynamic reductions is effective 
in any dimension: in particular, (3 -|- 1)- and (5 -|- l)-dimensional integrable examples were 
uncovered. We recall, see that the requirement of the commutativity of the flows Q is 
equivalent to the following restrictions on their characteristic speeds: 



(10) 



i 7^ j, dj = d/dfij, no summation! Once these conditions are met, the general solution of the 
system Q is given by the implicit 'generalized hodograph' formula |4()j . 

v\R) = x + X\R)t + fi\R)y, (11) 

i = l,...,n. Here u*(R) are the characteristic speeds of the general flow commuting with ©, 
that is, the general solution of the linear system 



djV^ 9jA* djfj.^ 



— A-J — A* fi^ — 



(12) 
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Substituting u(i?^, i?") into © and using @ one arrives at the equations 

{A + ^I'B + A^^) diu = 0, (13) 

(no summation) implying that both A* and /i* satisfy the dispersion relation 

D{fi, A) = det(^ + nB + \E) = 0. (14) 

Thus, the construction of solutions describing nonlinear interactions of n planar simple waves 
consists of two steps: 

(1) Reduce the initial system ® to a pair of commuting flows © by solving the equations (|in|). 
H13|l for u(R), A*(R), /i*(R) as functions of the Riemann invariants R^, These equations 
are highly overdetermined and do not possess solutions in general. However, for integrahle sys- 
tems these reductions depend, modulo reparametrizations f^{R^), on n arbitrary functions 
of a single argument. Once a particular reduction of the form Q is constructed, the second step 
is fairly straightforward: 

(2) Solve the linear system 1)12(1 for i;'(R) and determine R^, as functions of t,x,y from 
the implicit hodograph formula ([111) . This step provides some extra n arbitrary functions. 
Therefore, solutions arising within this scheme depend on 2n essential functions of a single 
argument. 

Simple waves. For n = 1 we have u = u{R), thus, hodographs of these solutions are curves. 
The scalar variable R = R^ solves a pair of first order PDEs 

Rt = ^{R) Rx, Ry = fJ-{R) Rx 

which, in the one-component situation, are automatically commuting. The hodograph formula 
PHI takes the form 

f{R) = X + X{R)t + n{R)y (15) 

where f{R) is arbitrary. This formula shows that, in coordinates t,x,y, the surfaces R = const 
are planes so that the solution u = u{R) is constant along a one-parameter family of planes. 
Therefore, it is singular along the developable surface which envelopes this one-parameter family. 
Solutions of this type, known as planar simple waves, exist for all multi-dimensional quasilinear 
systems and do not impose any restrictions on the matrices A and B. 

Double waves. For n = 2 we have u = u.{R^,R^) so that hodographs of these solutions (known 
as double waves, or nonlinear interactions of two planar simple waves) are surfaces. Double 
waves, as well as simple waves, belong to the class of solutions with a 'degenerate hodograph'. 
In the context of multi-dimensional gas dynamics, they were extensively investigated in |39j . 
Notice that each hodograph surface corresponds to infinitely many solutions, indeed, one needs 
to solve the system Q for R^,R^ to make a hodograph surface into a solution. This step 
does not change the hodograph surface, it just specifies the dependence on t, x, y. The general 
solution of Q is given by the implicit hodograph formula 

v\R) = x + X\K)t + ^x\K)y, t;2(R) =x + A2(R)t + /i2(R)y. (16) 

Setting R^ = const, R? = const, one obtains a two-parameter family of lines (or, using 
differential-geometric language, a line congruence) in the 3-space of independent variables t, x, y. 
The corresponding solution u = u{R^,R'^) is constant along the lines of this two-parameter fam- 
ily. Therefore, it is singular along the focal surfaces of the congruence. 
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One can show that any 2x2 system @ possesses infinitely many two-component reductions 
of the form Q parametrized by two arbitrary functions of a single argument. The corresponding 
double waves depend on four arbitrary functions. On the contrary, for multi-component systems 
(m > 3) the requirement of the existence of double waves imposes strong restrictions on the 
matrices A and B. In Sect. 3 we prove our main result 

Theorem 2 The Haantjes tensor of an arbitrary matrix ^ is zero if and only if 

(i) the system ^ possesses double waves parametrized by four arbitrary functions of a single 

argument; 

(a) the characteristic speeds A*, fj,^ of two- component reductions are not restricted by any alge- 
braic relations other than the dispersion relation D{fj,, X) = 0. That is, for any point uq in the 
hodograph space and any two points {X^,fi^) and on the dispersion curve at uq one can 

find a two-component reduction whose characteristic speeds at this point take the values (A^, A^) 
and {ji^ , ^'^) , respectively. 

The condition (ii) is crucial: for instance, the equations of two-dimensional gas dynamics 
possess potential double waves parametrized by four arbitrary functions, however, the system 
is not diagonalizable, see Sect. 5. The explanation of this phenomenon lies in the fact that 
the characteristic speeds of commuting flows defining two-component reductions come from one 
and the same branch of the dispersion curve (which is reducible), so that the condition (ii) is 
violated. 

Thus, we have an easy-to- verify necessary condition for the integrability of multi-component 
multi-dimensional systems of hydrodynamic type. It should be emphasized that the condition 
of diagonalizability is necessary, but not at all sufficient for the integrability. In Sect. 2 we 
construct an example of a diagonalizable semi-Hamiltonian 3x3 system in 2 -|- 1 dimensions 
which is not integrable for generic values of parameters (Example 2). 

Multiple waves. The requirement of the existence of nontrivial 3-component reductions im- 
poses further constraints on A and B. These prove to be very restrictive and imply the existence 
of n-component reductions for arbitrary n jl3j . This phenomenon is similar to the well-known 
three-soliton condition in the Hirota bilinear approach which, generically, implies the existence 
of n-soliton solutions and the integrability. For two-component systems © the full set of con- 
straints imposed on A and B by the requirement of existence of three-component reductions 
was obtained in |14j . These constraints imply, in particular, that all 2-component integrable 
systems of the form © are necessarily symmetrizable in the sense of |2f)) (that is, possess three 
conservation laws of hydrodynamic type). Moreover, they necessarily possess a scalar pseudopo- 
tential playing the role of a dispersionless 'Lax pair'. We expect that both results generalize to 
a multi-component situation. 

We emphasize that, although the method of hydrodynamic reductions provides infinitely 
many (implicit) solutions parametrised by arbitrarily many functions of a single argument, the 
question of solving the initial value problem for integrable systems © remains open. We believe, 
however, that solutions describing nonlinear interactions of planar simple waves are locally dense 
in the space of all solutions of © (see ^| for a discussion of this issue in the context of 
the dispersionless KP equation). A detailed investigation of their breakdown and singularity 
structure is important for the analysis of the general Cauchy problem for multi-dimensional 
quasilinear systems. We hope that the combination of the available existence results (for the 
Cauchy problem for semi-Hamiltonian systems in the dimension one) with the decomposition of 
an integrable multi-dimensional system into a collection of commuting semi-Hamiltonian flows 
would lead to a new understanding of the Cauchy problem in many dimensions. 
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2 Examples 

In this section we list some examples of diagonalizable semi-Hamiltonian 3x3 systems of hy- 
drodynamic type in 2 + 1 dimensions. 

Example 1. Let us consider the so-called generalized Benney system |43j . 

at + {av)x = 0, vt + vvx + = 0, Wy + = 0, 

which reduces to the shallow water equations in the limit y = —x, w = a. In matrix form, we 
have Ux + A{u)ut + B{u)uy = 0, here u = {a,v,wY or, explicitly. 
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= 0. 

y 

We have verified that the Haantjes tensor of the corresponding matrix Q is zero. Since the 
system possesses four conservation laws 

at + {av)x = 0, vt + {v^/2 + w)^ = 0, 

Wy + ax = 0, {av)t + {aw + av'^)x + {up' /2)y = 0, 

the semi-Hamiltonian property is automatically satisfied (we have verified it independently using 
the tensor criterion from the Appendix). 

Looking for hydrodynamic reductions of the generalized Benney system in the form a = 
a{R}, ...,i2"), h = b{R^, w = w{R^, ...,i2"), where the Riemann invariants satisfy ©, 

one obtains 

diw = -{X' + v)d^v, dia = n\X' + v)div, (17) 
along with the dispersion relation 

~~{x^ + vy 

The compatibility condition didjw = djdiW implies 



= -7T^^- (18) 



didnV = — rdiv + -— ^ -dnV, (19) 



while the commutativity condition ^\()\ results in 

A-'' + V 

= /. " d^v. (20) 
The substitution of (^0]) into H19() implies the system for v{R) and A*(i?), 

a A* = ~^^. djv, didjv = — ^,^-^ divdjv. (21) 
J A* - AJ ^ ' ^ (AJ - A*)2 * ^ ^ ^ 

One can verify that the remaining compatibility conditions didja = djdia are satisfied identically. 
For any solution A*, v of the system 1)2 ip one can reconstruct a, /i* by virtue of (fTTj) . ((TH|) . The 
system (|21|) is compatible, with the general solution depending on n arbitrary functions of a 
single argument (modulo reparametrizations i?* — > f^{R^)), thus manifesting the integrability of 
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the generalized Benney system. We also recall that the generalized Benney system arises as a 
consistency condition of the dispersionless Lax pair |43j . 

1 2 





7? 

^ uu 


^ uv 


TP 

^ uw 


^ uv 


TP 

^ vv 


TP 

^ vw 




TP 

^ vw 


TP 

^ WW 





': 








^ w 1 



2^^ ' i^.-v 
Two-component reductions are governed by the equations 

whose general solution depends, modulo reparametrizations B} f^{R^), — > f'^{E?'), on two 
arbitrary functions of a single variable. Two extra arbitrary functions come from the solution 
of the corresponding hydrodynamic type systems Since the characteristic speeds (A^, /i^) 

and (A^, fi^) are not restricted by any algebraic relations other then the dispersion relation H18() . 
both conditions of Theorem 2 are satisfied. 

Example 2. Let us consider a class of symmetric systems © for which A is a constant 
diagonal matrix and i3 is a Hessian: 



(22) 



Here c < b < a are constants and the potential F(u, v,w) is of the form 

F = v{-fu + 6w) + f{-fu + 6w) (23) 

where the constants 7 and 6 satisfy the constraint {b — a)6'^ + (6 — 0)7^ = and /(•) is an 
arbitrary function of a single argument (one can set 7 = 1/ \/b — c, 6 = 1/ \/a — b) . Potentials 
of this type arise in the classification of diagonalizable systems of the form (|22() - see Sect. 4. 
We have verified that the system (|22|) . (|23|) is diagonalizable and semi-Hamiltonian. In fact, 
the semi-Hamiltonian property follows from the conservative nature of ()22() . In spite of being 
diagonalizable and semi-Hamiltonian, these systems are not integrable in general (that is, do not 
possess n-wave solutions for n > 3). The integrability conditions impose one extra constraint 
on the potential (|23j) . 

f"" = i'^n2{f"r-f"f"")- 

c — a 

This ODE can be solved explicitly implying that, up to elementary changes of variables, 

/" = ^ cot (7^ + 5^1-), 6 = 7^^^, (24) 
y/e c — a 

provided 7^ 0; otherwise the system H22() is linear. Notice that we need /" (rather then / 
itself) to write down equations (|22|) : 

tV" 7 iSf" 

7 6 I I I • (25) 
j6f" 6 d^f" 
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For /" given by (pi]) this system possesses the scalar pseudopotential 

V't = + :^ log sin(7'u - pipy) + ^ log sm{6w + pipy) 

ipx = + ^log sin(7M - pipy) + ^ log sin(5u; + pip^ 



y)i 





(a 





0^ 









h 









^ 










where p = \^ (that is, the system (pSj) arises from the compatibility condition iptx = V'xt)- 
Thus, this pseudopotential can be viewed as the dispersionless analogue of the Lax pair. In the 
2x2 case, the existence of a scalar pseudopotential is the necessary and sufficient condition 
for the integrability of a multi-dimensional quasilinear system ^1]. This should be true in a 
multi-component situation as well, although the proof meets technical difficulties. 
Notice that the dispersion relation for the system ()25|). 

l^f 7 l^f 

7 (5 
7<5/" (5 /" 

factorises into a line and a conic, 

(A - 6m)[(A - a/.)(A - c/x) - (7^ + + /"(A - 6^))] = 0. 

Thus, the system ()25() provides an integrable example from the class discussed in jSj. 

Below we list some further examples of three-component integrable (2 + l)-dimensional sys- 
tems of the form Q which were constructed in P] using the classical i?-matrix approach (see also 
|29j ) . We have verified that all of them are diagonalizable and semi-Hamiltonian by directly com- 
puting the corresponding tensors © and H45() . In fact, the semi-Hamiltonian property follows 
from their conservative nature. For each of these examples we have calculated all conservation 
laws of hydrodynamic type and scalar pseudopotentials which play the role of dispersionless 
Lax pairs. We hope that these examples will be useful for a further research in the area of 
multi-dimensional systems of hydrodynamic type. 

Example 3. The 3-component system 
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-v/2 
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0^ 















II 




1/2 





0/ 







constitutes the 'second' commuting flow in the dispersionless KP hierarchy. It possesses four 
conservation laws of hydrodynamic type, 

ut = 2vx, 
vt = {2w-u'^/2)x, 

{2W - u'^/4:)t = Uy- {UV)X, 

{4uw + 2v'^ - u^/2)t = {u\ + {Svw - Sn^t;)^, 
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(i) 











and the scalar pseudopotential 

The system arises from the compatibiUty condition tpty = il)yt. 
Example 4. The 3-component system 

2-r \/^\ /lOOWn 

rv {l-r)u 2-r + u 

{l + r)w {l-r)u I \ w I ^ \ / \ u; 

r=const, possesses four conservation laws 

ut = Uy + {2- r)vx, 

^■u;i+^^ ={l — r)(^w^uj , 
{u' + = + ^((2 - r)w + ruv),, 

( l-2r \ / l-2r 2-r \ 

u; 1+r = 1^(1 — r)'Ui;ii; !+'■ + (1 + r)ii;i+'' J . 

and the scalar pseudopotential 

V't = (1 - r) I nV-x + -0^ 



2-r 



V'j/ = (1 - r) ( Wlpx + Vlpi ^ + Wtpx ^ + -01 



r+1 2-r 
r-1 I „/, 1-r 







(i) 











Example 5. The 3-component system 

{r-l)v {2-r)u \ / ^ \ / 1 \ / u 
(2 - r)u f + f 

1 + r I \ w I ^ \000/Vw; 

r=const, possesses four conservation laws 

wt = {1 + r)ux, 

( ''-1 \ / r-l \ / r-1 \ 

|^^2-r = |^l(2-r j + _ 1) 2-r ^; j ^ 

(i; + ^w;2)^ = (2-r)(t^t(;)^, 
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and the scalar pseudopotential 



■0i = (1 - r)u'ipx , 



r-1 
X 



3 Proof of Theorem 2 

Let us first clarify where the obstruction to the existence of double waves comes from. To 
construct double waves one needs to solve the commutativity equations (|lf 



along with the relations ()13|). 

{A + fi'B + X'E) diU = 0, 
here i,j = 1,2 and u = u(i?^,i?^). The last condition implies the dispersion relation 

det{A + fx'B + A^^) = D{n\ X') = 0. 

We assume that the matrix A + + XE has a simple spectrum for generic ^ and A satisfying 
the dispersion relation. Then relations ()13() give expressions for the derivatives of w*, s > 2, in 
terms of the derivatives of = u, 

diu' = A/r{fi\X^)diu, d2u' = M%fi'^,X'^)d2U, s = 2,...,m, (26) 

where M'^(/x,A) are rational functions of and A whose coefficients are determined by A and 
B. Applying the differentiation dj, j ^ i, to the dispersion relation Z?(/i*,A*) = 0, taking into 
account and the commutativity conditions (fTU)) . one obtains 

^ 1 = T2 TT = > > ^ > )^2n, 2 = Ti rn = Pif^ ,fJ' )diu (27) 

— fJ: A — A fJ, — /i A — A 

where, again, P is rational in its arguments. Finally, the consistency conditions for the equations 
(|26|l . did2U^ = d2diu^, imply m — 1 (a priory different) expressions for did2u: 



did2U = R\fi^,^j?,X^,X^)diud2U. 

Thus, the condition of the existence of double waves is 

i?"(/x\;ti2^A\A2) = i?'(/i\^^A\A2) mod Z)(/z\ A^) = ^(/i^, A^) = 0. 

Since coefficients of the rational expressions W , i?' depend on the first derivatives of the matrix 
elements of A and B, this constitutes a set of nontrivial first order constraints for A and B. 

An informal way to demonstrate the necessity part of Theorem 2 is the following: all of the 
above formulae possess the specialization = k+l^^ , k,l = const (this is actually a subtle point 
which requires a justification). Then the dispersion relation reduces to det{kE+A+iJL^[lE+B)) = 
0, while the relation (fT3|) takes the form {kE + A + fi^{lE + B)) diU = 0. Thus, diu and 82^ are 
the eigenvectors of the matrix (IE + B)~^{kE + A) corresponding to the eigenvalues /i^ and /i^. 
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Since the surface u{R^,R^) is tangential to both eigenvectors, the two-dimensional distribution 
spanned by them is automatically holonomic. Since this is true for any two eigenvectors of the 
matrix (IE + B)~^{kE + A)^ it is diagonalizable and, hence, has zero Haantjes tensor. A more 
rigorous proof of this result is given below. 

Theorem 2 The Haantjes tensor of an arbitrary matrix ^ is zero if and only if 

(i) the system |^ possesses double waves parametrized by four arbitrary functions of a single 

argument; 

(a) the characteristic speeds A% ^* of two-component reductions are not restricted by any alge- 
braic relations other than the dispersion relation D{^.,\) = 0. That is, for any point uq in the 
hodograph space and any two points {X^^fi^) and (A^,;U^) on the dispersion curve at uq one can 
find a two-component reduction whose characteristic speeds at this point take the values (A^, A^) 
and{^^,fj?), respectively. 

Proof: 

Let the system © possess sufficiently many double waves satisfying the conditions (i), (ii) 
of the Theorem. To prove the necessity part of Theorem 2 we need to show that the Haantjes 
tensor of any matrix Q equals zero. Let us fix uq in the hodograph space and consider a line 
X = k -\- Ifi in the (A, ^)-plane; here k, I = const are the same as in (jH)). The intersection of this 
line with the dispersion curve D{fi, A) = at uq consists of m points (X^,^^), {X^,fi^). Let 
us fix any two of them, say, (A^,^-*^) and (A^,^^). According to our assumptions there exists a 
double wave solution passing through uq such that the characteristic speeds of the corresponding 
two-component reduction at the point uq take the values (A^, A^) and (/x^,//^), respectively. We 
thus have 

{A + fi'B + X'E) diu = 0, 

i = 1,2. Setting = K -\- Lfj,^ (notice that at the point uq the values of K and L coincide with 
k and I, respectively) we obtain 

{KE + A + fi\LE + B)) diU = 0. 

Introducing V = -{LE + B)-^{KE + A) we arrive at VdiU = fi^diU or, using indices, 

V;^diu'' = fi'diu""; (28) 

(no summation over i). Substituting A* = X + L/x* into the commutativity conditions ()1U() we 
obtain the relations 

fi'djL + djK = 0, (29) 
^ / J) hj = ^i 2- Applying to (|28|) the operator dj, j / i, we obtain 

Vj^kdiu'^dju^ + V^didju"" = dju'diu"" + n'didju"". (30) 

It is important to emphasize that the derivatives of K and L in the left hand side cancel out 
by virtue of H29|) . Thus, both K and L behave like 'constants' in all tensor formulas below. 
Interchanging i and j in (|3()|) and subtracting the results we arrive at the expression for didjU^ 
in the form 

pi i r) j T/™ 

s-a^,- = ^^ii^a.u™ + ^^a-tx™ + — Ifidiu^^djuK (31) 
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Substituting this back into (|Hn|) we obtain a simple relation 



where is the Nijenhuis tensor of V [K and L can be regarded as 'constants' in the computation 
of A^). This can be rewritten in the invariant form 

djij^'diu + difi^dju = — -{—jj' — (32) 
which implies the following four relations: 



fi^ fiWj fi^ diU + {n^YdiH^ djU 







VN{Vdi 


u, dju) 






VN{diU 


, Vdju) 






N{V^^u, 


Vdju) 



For instance, the first relation can be obtained by applying the operator V'^ to (|H2|) and using 
VdiU = ^*9iU, etc. Taking a linear combination of the above relations we obtain 

V'^N{diU, dju) - VN{VdiU, dju) - VN{diU, Vdju) + N{VdiU, Vdju) = 0. 

The last formula can be rewritten in the form H{diU, dju) = where H is the Haantjes tensor, 
indeed, a coordinate-free form of Q is 

H{X, Y) = V^N{X, Y) - VN{VX, Y) - VN{X, VY) + N{VX, VY) 

where X, Y are arbitrary vector fields. Thus, the value of H on any pair of eigenvectors of the 
matrix V equals zero (the choice i = 1, j = 2 is not essential). Hence, the Haantjes tensor of 
the matrix V = —{IE + B)~^{kE + A) equals zero at the point uq, so that the Haantjes tensor 
of {kE + A)~^{IE + B) is zero, too. Since uq is arbitrary, this finishes the proof of the necessity 
part of Theorem 2. 

The above considerations can be readily inverted to establish the sufficiency part of the the- 
orem. To construct double waves one needs to solve the equations (|28|) . (|29|) . As demonstrated 
above, the consistency conditions for (|28() reduce to and ()32|) . respectively. Among the 
second order equations (|31|) only one is really essential, say, the equation for the first component 

of the vector u. The other are satisfied identically modulo H28|) and ()32() . It remains to point 
out that the vanishing of the Haantjes tensor is equivalent to the condition that the relation 
l|32|) is an identity: first of all, the vanishing of the Haantjes tensor implies that N(diU, dju) 
belongs to the span of diU and djU. Using the invariant definition of the Nijenhuis tensor 

N{X, Y) = [VX, VY] + [X, Y] - V[X, VY] - V[VX, Y] 

one obtains that the coefficients at 5jU and djU in both sides of (|.S2|) coincide identically. Thus, 
two-component reductions are governed by two first order relations H29() and one second order 
equation for u^. Up to reparametrizations — > f^{R^), — > f'^{R?') this leaves two arbitrary 
functions of one variable. Solving the corresponding hydrodynamic systems gives two extra 
arbitrary functions. Therefore, double waves depend on four arbitrary functions. Since K and 
L are non restricted by any algebraic relations, the condition (ii) is also satisfied. 
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4 A class of diagonalizable Godunov's systems 





(a 














h 









^ 





c) 





In this section we discuss a class of conservative 3x3 systems 

ut + ^(u)u^ + B{vi)viy = 

where u = {v} ^u^ ,v?y , ^ is a constant diagonal matrix and i? is a Hessian: 

Fii Fi2 -Fi3 
F21 F22 F23 
F31 F32 

here c < b < a are constants, Fij = d'^F/du^duK These systems belong to the class introduced 
in [201. In the 2x2 case the classification of integrable systems of this type was given in |14j . 
Here we concentrate on the diagonalizability aspect (which is automatically satisfied in the 2x2 
case). Calculating the Haantjes tensor W^j^ for the matrix {kE + A)~^{IE + B) and equating it 
to zero one arrives at an over-determined system of third order PDEs for the potential F. The 
simplest way to obtain these equations is the following. One first calculates the components Hl^ 
(two of the indices coincide) and writes each of them with common denominator. Numerators 
thereof are polynomials in k and / which are required to be identically zero. Setting successively 
k = —a, k = —b and k = — c in the expressions for these polynomials one arrives at the first set 
of relations, 

-^123 = 0) 



-^111 (^ - c) = -Fi22(c - a) + Fi33(a - b), 
F222(c -a) = Fii2(6 - c) + ^332(0 - b), 
F333 (a-b) = Fii3 (5 - c) + F223 (c - 



-^lll-?^23 — (-?^12-^13)l) -^222-^13 — (-?^21-?^23)2, -^333-^12 — (-?^31-?*32)3 

Fni = ((c - a){F?,)^ + {a - b){F?,)^)/A, 
F222 = {{b - c){Fl^)2 + (a - 6)(F|3)2)/A, 
i^333 = {{b - c){Fl), + (c - a)(F|3)3)/A, 



a). 

(33) 



where A = (5 — c)Fii + (c — a)F22 + (a — 6)^33. Before proceeding we point out that the system 
H33|) possesses two integrals, one quadratic and one fourth order. 



(34) 



/ = A2 + 4(c - a)(c - b)Fl^ + 4(6 - a) (6 - c)F1^ + 4(a - 5) (a - c)F^^, 

J={c- b)Ff^Fl + {b- a)Fl^Fl + (a - c)FI^FI + F^FiAA; 

both / and J are constant by virtue of (|33j) . The further analysis splits into two cases. 
Case 1. At least one of the expressions 

(a-c)Fl, + ia-b)F^„ [b - a)Fl + {b - c)FI, [c - a)Fl + [c - b)Fl (35) 

equals zero (these expressions appear as denominators when one solves the system (|33|) for the 
third derivatives of F). Let us suppose, for definiteness, that (6 — a)i^23 + {b — c)-F2i ~ 
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(notice that, since c < b < a, the other two possibiUties lead to complex- valued solutions). 
Differentiating this constraint by and and taking into account (|33j) i one obtains F112 = 
-P233 = 0. The substitution into ((221)3 implies F222 = 0. The further integration of the system 
(|33|) shows that, up to elementary transformations of variables, one has the following expression 
for F: 

F = u^{-fu^ + 6u^) + /(7u^ + 6u^); 

here the constants 7 and 6 satisfy the relation (6 — a)6'^ + (6 — 0)7^ = and / is an arbitrary 
function of the indicated argument. We refer to the Example 2 of Sect. 2 for a detailed discussion 
of this case. 

Case 2. All expressions (|35|) are nonzero. In this case one can solve the relations (|33)1 for the 
third derivatives of F. Calculating the remaining components -f^j^, ^ 7^ i 7^ ^1 of the Haantjes 
tensor for the matrix {kE + A)^^{IE + B) and equating them to zero one arrives at the three 
cases: 

2a. All third order derivatives of F equal to zero. This corresponds to linear systems. 

2b. The integral J equals zero. This case can be eliminated by the further consistency analysis. 

2c. The integrals / and J satisfy the relation 

/2 + 64(a-6)(&-c)(c-a)J = 0. (36) 

Notice that solutions constructed in Case 1 satisfy this constraint. We will show that there 
exists no other solutions in this class. First of all we point out that the general (real-valued) 
solution to the first four equations (|33|) (which are linear in F) is given by the formula 

F = g{z) + g{z) + p{a) + q{z + z + a) + s{a + z)+s{a + z)+T 

where 

12 12 13 

U . U _ U . U U U 

z = ~|~ 1 z = — X = — 

and T is an arbitrary quadratic form in v}-,v?^v? . Without any loss of generality one can assume 
T = 2e ^^_|. . Substituting this ansatz for F into H34|) and keeping in mind the relation (|36() one 
arrives at the following set of functional equations for the functions G = g" , P = p" , Q = 
q", S = s": 

{G{z) + G{z) + P{a) +Qiz + z + a) + S{a + z) + S{a + z) + 2e)^ - 
{G{z)-G{z)f -{P{a)-Q{z + z + a)f -{S{a + z)-S{a + z)f = I/A, 

{G{z) - G{z)f (P(a) -Q{z + z + a)f + 
(G(z) - G{z)f {S{a + z)- S{a + z)f + 
(P(a) -Q{z + z + a)f {S{a + z) - S{a + z)^ + (38) 
2 (G(z) - G{z)) {P{a) -Q{z + z + a)) {S{a + z) - S{a + z)) x 
(G(z) + G{z) + P(a) +Q{z + z + a) + S{a + z) + S{a + z) + 2e) = P/6A. 

These functional equations can be solved explicitly as follows. Imposing the constraint G{z) = 
G{z) in the second equation (|38() one obtains {P{(y) — Q{z + z + a))^ [S{a + z) — S{a + z)) = 
I^/M, so that 

(P(a) -Q{z + z + a)) {S{a + z) - S{a + z)) = 1/8 (39) 
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(the case of the opposite sign is considered in a similar way). Setting G{z) = G{z) in the first 
equation (|37|) we have 

{G{z) + G{z) + P{a) + Q{z + z + a) + S{a + z) + S{a + z) + 2ef = 
{P{a) -Q{z + z + a)f + {S{a + z) - S{a + z)f + //4 = 
(P(a) - Q{z + z + a) + S{a + z) - S{a + z))^ 

by virtue of Thus, 

G{z) + G{z) + P{a) + Q(z + z + a) + S(a + z) + S{a + z) + 2e = 
-P{a) + Q{z + z + a) - S{a + z) + S{a + z) 

(the case of the opposite sign is considered in a similar way) so that G{z)+P{a)+S{a+z)+e = 0. 
This simple functional equation implies that G,G, P, S, S are linear functions: 

G{z) = cz + fi, G{z) = cz + fl, P{a) = ca + z^, 
S{a + z) = —c{a + z) + r], S{a + z) = — c(a + z) + f/. 

However, the substitution of these expressions into (|S7j) readily implies c = 0. Thus, G, S and 
P are constants. This brings us back to the Case 1 considered above. 

5 Double waves for two-dimensional gas dynamics 

The equations of two-dimensional isentropic gas dynamics are of the form 

Pt + ipu)x + {pv)y = 0, Ut + UUx + VUy + Px/P = 0, Vt+ UVx + VVy + Py/p = 0, 

where p = p{p) is the equation of state. In matrix form, one has + Aux + Buy = where 
u = (yO, u, vY and 

/ u p \ / V p \ 

A = c^/p u \ , B = I I'D ; 
\ n / \ cV/O V / 

here = p'{p) is the sound speed. We have verified that, although the Haantjes tensor of 
an arbitrary matrix from the linear pencil A + kB equals zero, this is not the case for the 
general family Q. In particular, the Haantjes tensor of the matrix A~^B does not vanish. 
Nevertheless, equations of gas dynamics possess double waves of special type, namely, potential 
double waves. These solutions have been extensively investigated in [23], see also references 
therein. We discuss them in a different setting below. Let us first recall that the dispersion 
relation det(£' + \A + pB) = factorises into a line and a conic, 

(1 + Au + pv){{l + An + pvf - c2(A2 + p'^)) = 0, 

the conical branch corresponding to sound waves. As demonstrated in 39 , there exist potential 
flows describing nonlinear interaction of two sound waves which are locally parametrized by 
four arbitrary functions of a single argument. To construct these solutions we use the ansatz 
p = p{R^,R^), u = u{R^,B?), V = v{R^, R?) where the Riemann invariants R^^R? solve a pair 
of diagonal systems 

i?; = V(R)i?^ < = ^'(R)i?j, (40) 
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i = 1,2. The substitution implies the relations 

(1 + X'u + fi'v)diu + ^X'dip = 0, 

(1 + X'u + p'v)div + fp'dip = 0, 

(1 + X^u + p^v)dip + p{X^diU + p^div) = 0, 

i = 1,2, di = d/dEJ'. Expressing diU and diV from the first two equations and substituting them 
into the third, one arrives at the dispersion relation for sound waves, 

(1 + x'u + p'vf - c\{X'f + {p'f) = 0. 

Parametrising A* and /x* in the form A* = cos 99*, /U* = s* sin ip^ , we obtain 1 + X^u + p^v = cs^ 
so that = l/(c — ncos(/j* — vsvuip"^). Thus, 

A^ = p^ = (41) 

c — u cos if^ — V sm c — u cos cp^ — v sm (p^ 

The equations for u and v take the form 

9iu + ^ cos <p^dip = 0, 92"" + ^ cos 'p'^d2P = 0, 

(42) 

c^iz; + ^ sin ip^dip = 0, 92ii + ^ sin (f'^d2P = 0. 

Notice that since p^diU = X^diV one has Uy = Vx- Thus, solutions describing nonlinear in- 
teraction of two sound waves are necessarily potential. Writing out the commutativity con- 
ditions djX^/{X^ — A*) = djp^/{p^ — p^), i,j = 1,2, i / j, and the consistency conditions 
did2U = d2diu, did2V = d2div, one arrives at the following system for ip^, ip'^, p: 



02(^1 = cot ^ (4 + i cos ((^1 - (^2) j d^p^ 
ai^2 ^ + i cos (VPI - ^^)) dip, 

did2P = (- COS (^^ - ^^) + ^(3 + COS 2(^1 - c,^))) 



(43) 



Up to reparametrizations f^{R^), the general solution of this system depends on two 

arbitrary functions of a single argument. For any solution 93* (i?^, i?^), p{R^,R^) one recon- 
structs the hodograph surface of the corresponding double wave by solving the equations 
for u{R^,R^) and v{R^,R?) (which are automatically consistent). Each of these surfaces can be 
made into a solution of the gas dynamics equations by solving the equations (|in|) for R^{t, x, y) 
and R'^{t, x, y); here the characteristic speeds A*, /n* are given by (|1T|) . Notice that since the gen- 
eral solution of (|4()|) depends on two arbitrary functions of a single argument, each hodograph 
surface corresponds to infinitely many solutions. 

Remark. Any hodograph surface defined parametrically in the form u{R^,R?), v{R^ , R?), 
p{R^,R?) can be parametrized explicitly as p = p{u,v). Using the relations 1)42(1 . (|43|) one can 
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show that the function p{u, v) satisfies the second order PDE 

h i^) {pl + pI)-^ Uipl + pI) + 2^ - {pi + pD/p) = 

(44) 

~7l{Puu ^ Pvv) PvPuu ~l~ '^PuPvPuv PuPvv) 

here = p' where p{p) is the equation of state and ' = d/dp. For a polytropic gas, p = p'^ J'^-, 
the equation simphfies to 

(3 - l)p^~''^{p'i + pI) - 2p^~'^"' = p^~'^{puu + Pvv) - pIpuu + 2puPvPuv - pIpvv 
which takes the form 

(1 + iPuu){^l - (7 - - 2(puipv^uv + (1 + ^vv){^l - (7 - 1)95) = 

after the substitution ip = p'^~^ /{'J — 1)- In this form the equation for double waves appears in 

m, §8. 

It is unhkely that the system (|43|1 can be integrated in a closed form, even for Chaplygin's 
equation of state p{p) = b — a'^/p, a,b = const, in which case it simplifies to 

a^^i = sin (^1 - y,^), dw^ = ^ sin ((^2 _ ^1)^ 

9^^2/0 = 2^i^sin2£i^. 

However, some particular solutions can be constructed. Let us consider the polytropic case, 
p{p) = p'^ , so that ^ = One can verify that 

(f = R — R , (f = R — R + —, p = exp t-i^ ' 
solve the system (|43() . The corresponding equations imply 

u = ^(cos(iii - - sin(i?i - i?2))exp-(^'+^'), 

V = J^(cos(i?i - + sin(i?i - exp-(^'+«') . 

Excluding R^ and R^ one arrives at the explicit hodograph surface, ^^^^^(u^ + v'^)p^~^ = 1, or, 
equivalently, (1 — 'y)'^{u'^ + u^) = 20^. To make this surface into a solution one has to solve the 
equations (|4flj) where the characteristic speeds are given explicitly by (|41|) : 

= ^ cos(i?i - exp^^+«\ A2 = sin(i?i - exp«^+«^ 

/.I = ^ sin(i?i - exp«^+^^ ^2 = l_i cos(i?i - ii^) exp«^+«^ 

The general solution R^{t, x,y), R?{t,x,y) of the corresponding system 1)40(1 is given by the 
implicit formula 

^^(R) = t + A^(R)x + Ai^(R)y, t-^CR) = t + \^{K)x + ;u2(R)y 
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where v^,v'^ are the characteristic speeds of the general flow commuting with ()41|1. They satisfy 
the linear equations 

Potential double waves of this type and their applications to problems of gas dynamics were 
extensively investigated in the monograph [HHI • The aim of this section was to demonstrate how 
solutions with a degenerate hodograph fit into the general scheme of the method of hydrodynamic 
reductions proposed in ^Hl- In particular, our derivation of equations governing potential double 
waves is, in our opinion, much more transparent and straightforward since the use of Riemann 
invariants R^,B? considerably simplifies the calculation of compatibility conditions. On the 
contrary, the approach of I3Q, leads to quite cumbersome calculations based on Cartan's theory 
of over-determined systems. 

6 Appendix: the semi-Hamiltonian property 

For diagonalizable one-dimensional systems there exists a tensor object responsible for the 
semi-Hamiltonian property. We emphasize that the method does not require the actual trans- 
formation to the diagonal form ((21): it is sufficient to verify that the corresponding Haantjes 
tensor (jSJ is zero. After that one computes the (1, 3)-tensors M and 

MLi = mpvP,^ + K^vlN^^ - N;^Nf,v] - N;^NI.vI - Ni^Nf^v] - iV|,iV^>,' 

and 

+bid^Nl - b]dur>Ni^ + Nl,dur>b<^^ - iv;pa„.6f + iV;/„,6P; 

here h = v"^, that is, 6*- = VpV^. Using M and K one defines the (l,3)-tensor Q as 
Finally, one introduces the tensor P, 

PL, = 4QU< + ^'vQU^' - ^HQl, - QlvA^V (45) 

Theorem 3 '^33] A diagonalizable hydrodynamic type system with the matrix Vj{u) with 
mutually distinct eigenvalues is semi-Hamiltonian if and only if the corresponding tensor P is 
identically zero. 

Note that these objects can be obtained using computer algebra. 



19 



Acknowledgements 



We thank M.V. Pavlov and S.P. Tsarev for numerous discussions. We also thank the referees 
for a constructive criticism and useful suggestions. 

References 

M. Blaszak, B.M. Szablikowski, Classical R-matrix theory of dispersionless systems: II. 
(2+l)-dimension theory, J. Phys. A 35 (2002) 10345-10364. 

M. Burnat, The method of Riemann invariants for multi-dimensional nonelliptic system, 
Bull. Acad. Polon. Sci. Sr. Sci. Tech. 17 (1969) 1019-1026. 

M. Burnat, The method of Riemann invariants and its applications to the theory of plas- 
ticity. I, II. Arch. Mech. (Arch. Mech. Stos.) 23 (1971), 817-838; ibid. 24 (1972), 3-26. 

M. Burnat, The method of characteristics and Riemann's invariants for multidimensional 
hyperbohc systems, Sibirsk. Mat. Z. 11 (1970) 279-309. 

S. Canic and B. L. Keyfitz, A useful class of two-dimensional conservation laws. Proceedings 
of ICIAM 95: Supplement 2: Applied Analysis, Mathematical Research, 87, (K. Kirch- 
gassner, O. Mahrenholtz and R. Mennicken, editors) Akademie Verlag, Berlin, ZAMM 
(1996) 133-136. 

C. Dafermos, Hyperbolic conservation laws in continuum physics. Springer- Verlag, 2000. 

L. Dinu, Some remarks concerning the Riemann invariance, Burnat-Peradzyhski and Martin 
approaches. Rev. Roumaine Math. Pures Appl. 35, N3 (1990) 203-234. 

B.A. Dubrovin and S.P. Novikov, Hydrodynamics of weakly deformed soliton lattices: dif- 
ferential geometry and Hamiltonian theory, Russian Math. Surveys 44 , N6 (1989) 35-124. 

B.A. Dubrovin, Geometry of 2D topological field theories, Lect. Notes in Math. 1620, 
Springer- Verlag (1996) 120-348. 

E.V. Ferapontov, D. A. Korotkin and V.A. Shramchenko, Boyer-Finley equation and sys- 
tems of hydrodynamic type. Class. Quantum Grav. 19 (2002) L205-L210. 

E.V. Ferapontov and M.V. Pavlov, Hydrodynamic reductions of the heavenly equation. 
Class. Quantum Grav. 20 (2003) 2429-2441. 

E.V. Ferapontov and S.P. Tsarev, Systems of hydrodynamic type that arise in gas chro- 
matography. Riemann invariants and exact solutions. Mat. Model. 3 , N2 (1991) 82-91. 

E.V. Ferapontov and K.R. Khusnutdinova, On integrability of (2+l)-dimensional quasilin- 



ear systems. Comm. Math. Phys. 248 (2004) 187-206; arXiv:nhn.SI/03Q5044 



E.V. Ferapontov and K.R. Khusnutdinova, The characterization of two-component (2+1)- 
dimensional integrable systems of hydrodynamic type, J. Phys. A: Math. Gen. 37, N8 
(2004) 2949 - 2963; jarXiv:nlin.SI/031002l| 



20 



[15] E.V. Ferapontov and K.R. Khusnutdinova, Hydrodynamic reductions of multi-dimensional 
dispersionless PDEs: the test for integrability, J. Math. Phys. 45, N6 (2004) 2365-2377; 
|arXiv:nhn.SI/0312015, 

[16] E.V. Ferapontov and D.G. Marshall, Differential-geometric approach to the integrability of 
hydrodynamic chains: the Haantjes tensor, arXi v:nlin.SI/ 0505013( 

[17] J. Gibbons and Y. Kodama, A method for solving the dispersionless KP hierarchy and its 
exact solutions. II. Phys. Lett. A 135, N3 (1989) 167-170. 

[18] J. Gibbons and S. P. Tsarev, Reductions of the Benney equations, Phys. Lett. A 211 (1996) 
19-24. 

[19] J. Gibbons and S. P. Tsarev, Conformal maps and reductions of the Benney equations, 
Phys. Lett. A 258 (1999) 263-271. 

[20] S.K. Godunov, An interesting class of quasi-linear systems, Dokl. Akad. Nauk SSSR 139 
(1961) 521-523. 

[21] A. Grundland and R. Zelazny, Simple waves in quasilinear hyperbolic systems. I, II. Rie- 
mann invariants for the problem of simple wave interactions. J. Math. Phys. 24, N9 (1983) 
2305-2328. 

[22] F. Guil, M. Manas and L. Martinez Alonso, On the Whitham Hierarchies: Reductions and 
Hodograph Solutions, J. Phys. A 36 (2003) 4047-4062. 

[23] J. Haantjes, On X^-forming sets of eigenvectors, Indagationes Mathematicae 17 (1955) 
158-162. 

[24] I.M. Krichever, The averaging method for two-dimensional "integrable" equations, Funct. 
Anal. Appl. 22, N3 (1988) 200-213. 

[25] I.M. Krichever, Spectral theory of two-dimensional periodic operators and its applications, 
Russian Math. Surveys 44, N2 (1989) 145-225. 

[26] I.M. Krichever, M. Mineev-Weinstein, P. Wiegmann and A. Zabrodin, Laplacian growth 
and Whitham equations of sohton theory, Phys. D 198, no. 1-2 (2004) 1-28. 

[27] A. Majda, Compressible fluid flows and systems of conservation laws in several space vari- 
ables, Appl. Math. Sci., Springer- Verlag, NY, 53 (1984) 159 pp. 

[28] M. Manas, L. Martinez Alonso and E. Medina, Reductions and hodograph solutions of the 
dispersionless KP hierarchy, J. Phys. A: Math. Gen. 35 (2002) 401-417. 

[29] M. Manas, 5'-functions, reductions and hodograph solutions of the r-th dispersionless mod- 
ified KP and Dym hierarchies, J. Phys. A 37, no. 46 (2004) 11191-11221. 

[30] L. Martinez Alonso and A.B. Shabat, Hydrodynamic reductions and solutions of a universal 
hierarchy, Teoret. Mat. Fiz. 140, no. 2 (2004) 216-229. 

[31] M.V. Pavlov, Integrable hydrodynamic chains, J. Math. Phys. 44 (2003) 4134-4156. 



21 



[32] M.V. Pavlov, Classification of the integrable Egorov hydrodynamic chains, Theor. and 
Math. Phys. 138, no 1 (2004) 55-70. 

[33] M.V. Pavlov, S.I. Svinolupov and R.A. Sharipov, An invariant criterion for hydrodynamic 
integrability, Funktsional. Anal, i Prilozhen. 30 (1996) 18-29; translation in Funct. Anal. 
Appl. 30 (1996) 15-22. 

[34] Z. Peradzynski, Riemann invariants for the nonplanar /c-waves. Bull. Acad. Polon. Sci. Sr. 
Sci. Tech. 19 (1971) 717-724. 

[35] Z. Peradzynski, Nonlinear plane fc-waves and Riemann invariants, Bull. Acad. Polon. Sci. 
Sr. Sci. Tech. 19 (1971) 625-632. 

[36] B. Sevennec, Geometrie des systemes hypcrboliqucs dc lois de conservation, Memoire (nou- 
velle serie) N56, Supplement au Bulletin de la Societe Mathematique de Prance 122 (1994) 
1-125. 

[37] D. Serre, Systems of conservation laws. 1. Hyperbolicity, entropies, shock waves, Cambridge 
University Press (1999) 263 pp. 

[38] D. Serre, Systems of conservation laws. 2. Geometric structures, oscillations, and initial- 
boundary value problems, Cambridge University Press (2000) 269 pp. 

[39] A.F. Sidorov, V.P. Shapeev and N.N. Yanenko, The method of differential constraints and 
its applications in gas dynamics, "Nauka", Novosibirsk (1984) 272 pp. 

[40] S.P. Tsarev, Geometry of hamiltonian systems of hydrodynamic type. Generalized hodo- 
graph method, Izvestija AN USSR Math. 54, N5 (1990) 1048-1068. 

[41] S.P. Tsarev, Diagonalizable averaging of the generalized KdV equation, Uspekhi Matem. 
Nauk. 46, N6 (1991) 194. 

[42] L. Yu, Reductions of dispersionless integrable hierarchies, PhD Thesis, Imperial College, 
London, 2001. 

[43] E.V. Zakharov, Dispersionless limit of integrable systems in 2 -|- 1 dimensions, in Singular 
Limits of Dispersive Waves, Ed. N.M. Ercolani et al.. Plenum Press, NY, (1994) 165-174. 



22 



